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ABSTRACT 

One  Dimensional,  Finite  Time  Release  Spill  Model  With  Decay  Function 

Program  code  name:  SPILL 


Developed  by:  Wm.  C.  Allison  Distributed  by:  Z.  Novak 

Waterloo  Co-op  Student       River  Systems  Unit 
River  Systems  Unit  Water  Resources  Branch 

Water  Resources  Branch       Ministry  of  the  Environment 
Ministry  of  the  Environment   135  St.  Clair  Avenue  West 

Toronto,  Ontario 


This  is  the  original  micro-computer  version  of  the  program  completed 
on  March  30,  1985.  The  spill  concentration  equations  are  based  on  a 
"mainframe"  FORTRAN  version  of  the  spill  model. 

The  source  language  of  this  program  is  "BASICA".  The  executionable 
program  is  in  "COMPILED  BASIC-VERSION  1.0".  An  IBM  PC  or  compatible 
clone  is  required  to  execute  this  program.  The  micro-computer  must 
contain  at  least  256  K  of  memory. 

The  program  has  been  extensively  checked  for  coding  and  logic  errors 
by  the  developer  and  River  Systems  Staff.  The  model  has  been 
compared  to  a  documented  spill  in  the  Speed  River  with  good 
results.  The  details  of  these  tests  are  outlined  in  the  Engineering 
and  Software  documentation  Section  7.5. 

The  model  is  based  upon  the  "Fickian"  concept  of  spill  dispersion. 
This  version  allows  for  variations  in  the  mass  flux  of  the  spill 
over  a  finite  period  of  time.  The  model  assumes  that  the  spill 
originates  from  a  planar  source.  Only  longitudinal  dispersion  from 
the  spill  cloud  is  considered  and  this  is  accounted  for  by  using  an 
experimentally  or  analytically  derived  dispersion  coefficient.  An 
exponential  decay  function  may  be  utilized  if  the  substance  is 
non-conservative  and  a  decay  rate  is  known.  Refer  to  the 
Engineering  and  Software  Documentation  and  the  Users  Manual  for 
detailed  information  concerning  these  assumptions  and  the 
engineering  theory  the  model  is  based  upon. 


2.0  DESCRIPTION  OF  ENGINEERING  PROCEDURE 

The  model  calculates  downstream  concentrations  of  a  substance  that 
has  been  spilled  into  a  uniformly  flowing  watercourse.  Two  options 
are  available  in  the  model;  discrete  concentrations  at  specific 
times  and  distances  downstream  or  a  concentration  profile  over  time 
at  a  specified  distance  downstream. 

This  model  uses  modified  "Fickian"  spill  dispersion  theory  to 
calculate  downstream  concentrations  of  a  spill.  Changing  mass 
flux  rates  of  a  spill  are  accounted  for  by  superimposing  the  time  of 
each  flux  rate.  An  exponential  decay  rate  may  be  applied  to  the 
spill  if  appl icable. 

The  "Fickian"  theory  of  spill  dispersion,  which  this  model  is  based 
upon,  contains  many  assumptions.  The  validity  and  accuracy  of  much 
of  the  theory  has  been  disputed  (Section  3.4),  However,  this  model 
is  ideally  suited  as  a  preliminary  or  emergency  planning  tool  when  a 
fast  and  relatively  accurate  solution,  that  requires  a  minimum  of 
input  parameters,  is  desired.  A  potential  use  of  this  program  is  to 
profile  over  time  the  concentration  of  a  spill  passing  a  water 
intake  or  recreation  area.  Authorities  could  plan,  using  sufficient 
safety  factors,  when  to  discontinue  use  of  the  water  intake  or 
recreation  areas. 

Authorities  could  also  use  the  model  to  locate  a  spill  cloud  for 
possible  clean-up  action.  Downstream  concentrations  after  the 
effects  of  dilution,  dispersion  and  decay  (if  applicable)  may  also 
be  calculated  to  enable  the  user  to  judge  the  severity  of  a  spill. 
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3.0  ENGINEERING  PROCEDURE 


3.1  Mathematical  Model 


This  model  is  based  upon  a  modified  version  of  the  "Fickian"  spill 
dispersion  model.  True  "Fickian"  spill  dispersion  theory  requires 
an  assumption  of  instantaneous  release  of  the  spilled  substance.  In 
reality,  spills  often  last  long  periods  of  time  and  have  changing 
mass  flux  rates.  This  model  accounts  for  these  eventualities  by 
superimposing  the  release  time  of  each  spill  flux.  Intermittant 
spill  stoppages  can  be  accounted  for  by  using  a  mass  flux  rate  of 
zero. 

The  optional  exponential  decay  rate  of  the  spilled  substance  is 
applied  after  each  downstream  concentration  is  calculated.  The 
decay  rate  equation  is  of  the  form: 

CN  =  CO  *  EXP  (-1*DECAY*T) 

CN New  Concentration  Rate  After  Application  of 

Exponential  Decay  Rate. 

CO: Original  Concentration. 

EXP Exponential  Power  =  2.718**  Decay 

Decay...  Decay  Rate 

T Relative  Time  From  the  Start  of  the  Spill  to  Each 

Concentration  Calculation. 

Ideally,  prior  to  a  spill,  a  dispersion  coefficient  would  have  been 
calculated  from  experimental  data  for  the  section  of  waterway 
downstream  from  a  potential  spill  source.  However,  this  is  not 
always  possible  as  most  spills  are  unpredictable.  An  analytically 
derived  dispersion  coefficient  that  best  matches  the  physical 
properties  of  the  waterway  must  be  estimated.  Section  7.4  of  the 
Engineering  and  Software  Documentation  and  Section  3.0  of  the  Users 
Manual  deals  specifically  with  calculating  dispersion  coefficients 
and  recommends  a  suitable  range  of  numbers  for  various  physical 
situations. 
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3.2  Capabilities  of  the  Engineering  Procedure 

The  preliminary  input  parameters  required  are;  name  of  substance 
spilled,  optional  decay  rate  application,  name  of  the  watercourse 
the  spill  occurred  in,  the  name  of  the  source  of  the  spill  and  the 
start  and  end  times  and  dates  of  the  spill.  This  information  is 
necessary  in  order  to  uniquely  identify  each  simulated  spill  for 
future  reference. 

Physical  characteristics  of  a  watercourse  such  as;  width,  depth, 
velocity  and  a  dispersion  coefficient  are   required.   It  must  be 
stressed  that  these  input  parameters  are  average  values  reflecting 
the  topography  and  physical  characteristics  of  the  simulated  reach. 

The  spill  history  must  be  entered  into  the  model.  A  simple  steady 
state  spill  requires  only  a  flux  and  a  time  of  release.  A  spill 
with  a  changing  mass  flux  rate  must  be  divided  into  portions  such 
that  each  mass  flux  is  approximately  steady  state  for  each  time  of 
release  increment. 

Although  safety  factors  are   not  incorporated  into  the  program,  the 
user  should  always  utilize  conservative  values  for  input 
parameters.  Similarly,  water  quality  objectives  should  be 
conservative  to  insure  a  suitable  safety  factor. 

The  model  produces  downstream  concentrations  of  the  substance  after 
it  has  been  subject  to  the  effects  of  dilution,  dispersion  and  decay 
(if  applicable).  The  user  has  the  choice  of: 

1)  one  or  more  discrete  concentrations  at  specified  times  and 
distances  downstream. 

2)  a  concentration  profile  over  time  at  a  specified  distance 
downstream  with  a  specified  time  increment  between 
concentration  calculations. 
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A  tabular  output  is  produced  that  displays;  identifying  information, 
concentration,  time  of  travel  and  the  distance  downstream  to  the 
concentration  calculation  (fig.  1). 

This  table  appears  on  the  screen  and  is  available  for  hardcopy. 
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Figure  1 
Sample  Output 


Discrete  Concentrations 


*Htt ********************** ****^^FILL    PROP  ILE********************* ********* 
************************************************************************** 


CN 


SUBSTANCE  SPILLED  WAS:  CYANIDE 

DECAY    RATE    IS    .OOiiiOOl    S^-I 

IN  THE  HUMBER  RIVER 

THE  SOURCE  OF  THE  SPILL  WAS:  MACDONALDS 

TIME  OF  SPILL:  FROM  2:30  P.M.  JAN.  11  1 


01  M6/L 


NUMBER 


1 

2 
3 
4 
S 


TIME  FROM  START 
OF  SPILL  (hrs) 

20.000 
10.000 
40. 000 
30 .  i!>00 
15. 000 


TO  6:30  P, 

,M. 

JAN 

11  1986 

CONCENTRATION 

DISTANCE 

mg/1 

m 

0. 121D-01 

1 0':>oo 

0.3220+00 

4500 

0.  lOBD+00 

11000 

0.377D+00 

12000 

0.6430+00 

6000 

Concentration  Profile 


**************************************** *****************************%**** 

*************************** ***SF ILL    PROF  ILE******************** ********** 
t*********t*************** ******************************************  ****** 

SUBSTANCE  SPILLED  WAS:  CYANIDE    CN    .01  MS/L 

DECAY  RATE  IS  .000001  S — 1 

IN  THE  HUMBER  RIVER 

THE  SOURCE  OF  THE  SPILL  WAS:  MACDONALDS 

TIME  OF  SPILL:  FROM  2:30  P.M.  JAN.  11  1986  TO  6:30  P.M.  JAN  11  19B6 

THIS  PROFILE  IS  AT  A  DISTANCE  OF   10000   meters  DOWNSTREAM  FROM  THE  SPILL 


NUMBER 

TIME  FROM  START 

OF 

SPILL  (hrs) 

1 

17.778 

2 

19.778 

3 

21.778 

4 

23.778 

S 

25.778 

6 

27.778 

7 

29.778 

8 

31.778 

9 

33.778 

10 

35.778 

n 

37.778 

12 

39.778 

13 

41.778 

14 

43. 778 

IS 

45.778 

16 

47.778 

CONCENTRATION 
mg/I 

0 ,  '!'Oi!>D+*!'0 
0.9270-02 
0.686D-01 
0.259D+00 
0.5710+00 
O. 8200+00 
0. 8290+00 
0. 6260+00 
0.374D+00 
0. 183D+00 
0. 7520-01 
0.269D-01 
O. 848D-02 
0.2760-02 
0. 7040-O3 
O.OOOD+00 
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3.3  Limitations  of  the  Engineering  Procedures 

This  model  is  applicable  to  most  types  of  watercourse  that  flow  in  a 
uniform  direction.  Allowances  can  be  incorporated  into  the 
dispersion  coefficient  to  allow  for  non-typical  characteristics  such 
as  stilling  basins,  rapids,  or  extremely  winding  rivers. 

Variances  between  laminar  and  turbulant  flow  can  be  allowed  for  by 
weighting  the  dispersion  coefficient  to  reflect  changing  conditions 
in  the  reach. 

The  physical  stream  characteristics  must  be  positive  real  numbers, 
greater  than  zero.  The  spill  input  characteristics  must  be  positive 
real  numbers,  however,  a  mass  flux  of  zero  grams  per  second  is 
allowable  in  order  to  account  for  intermittent  stoppage  of  the  spill 
flow.  The  time  increment  between  calculations  is  limited  in  size  to 
the  magnitude  of  the  user's  patience  (smaller  increments  require 
longer  time  for  the  model  to  execute)  or  a  maximum  of  200 
concentration  calculations  is  allowed. 

The  decay  rate  (if  applicable)  must  be  a  positive  real  number  or 
zero.  Negative  numbers  intended  to  reflect  a  growth  rate  are  NOT 
al lowed. 

The  numeric  magnitude  of  some  equations  may  exceed  the  capacity  of 
the  microcomputer.  These  instances  of  numeric  overflow  usually 
occur  when  the  concentration  is  calculated  close  to  the  source  (i.e. 
under  5  metres)  and  the  velocity  is  very  fast  (i.e.  greater  than  20 
m/s).  There  is  a  built  in  facility  to  circumvent  this  possibility, 
however,  it  is  not  foolproof  and  overflow  may  occur  (sect. 5. 2).   If 
overflow  occurs  the  user  should  verify  the  accuracy  of  the  input 
data.  The  reliability  of  the  model  cannot  be  guaranteed  when 
extreme  or  non-typical  input  data  are  used.  The  model  is  suited  for 
typical  Ontario  Rivers  (Refer  to  section  7.5). 
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3.4  Engineering  Assumptions 


The  "Fickian"  spill  dispersion  theory  contains  three  simplifying 

3 

assumptions: 


1)  The  spill  originates  from  a  plane  source; 

2)  Only  longitudinal  dispersion  occurs; 

3)  That  concentration  versus  time  is  normally  distributed. 

The  first  and  second  assumptions  are  valid  if  concentrations  are 
calculated  below  a  point  downstream  of  complete  uniform  natural 
lateral  dispersion.  If  concentrations  are  calculated  before  the 
point  of  total  uniform  lateral  dispersion  the  model  will  produce 
concentrations  lower  than  what  they  actually  are. 

There  are  many  conflicting  opinions  on  how  to  estimate  the  point  of 
complete  lateral  mixing.  It  is  best  to  use  common  sense  when 
judging  the  validity  of  results  calculated  near  the  source  of  the 
spill.  There  have  been  a  number  of  mathematical  equations  proposed 
to  predict  this  point,  one  equation  is  presented  in  Appendix  1.0. 

Field  studies  have  shown  that  the  concentration  versus  time  profile 
for  a  spill  in  a  natural  watercourse  is  skewed.  Therefore,  the 
third  assumption  of  an  ideal  normally  distributed  profile  is  not 
valid.  The  user  should  apply  sufficient  safety  factors  to  account 
for  a  possible  offset  of  the  concentration  profile  calculated  by  the 
model . 
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4.0  DESCRIPTION  OF  THE  COMPUTATIONAL  PROCEDURE 

This  program  relies  entirely  on  input  data  to  simulate  a  spill  in  a 
waterway.  The  initial  information  required  is;  name  of  substance 
spilled,  option  of  decay  function,  magnitude  of  the  decay  function, 
name  of  the  waterway,  name  of  the  source,  start  and  end  times  and 
dates.  This  information  is  necessary  to  uniquely  identify  and  label 
each  spill.  The  waterway  characteristics  that  are  required  are; 
width,  depth,  velocity  and  a  dispersion  coefficient.  The  spill 
release  times  and  mass  fluxes  must  also  be  entered. 

The  user  has  two  output  options;  discrete  concentrations  at  a  series 
of  times  and  distances  downstream  or  a  concentration  profile  over 
time  at  a  point  downstream.  The  concentrations  for  both  options  are 
calculated  the  same  subroutines. 

The  concentration  calculation  subroutine  checks  to  insure  that  the 
spill  has  ended  at  the  time  the  concentration  is  calculated.  If  the 
spill  is  still  occurring,  the  spill  time  and  mass  flux  is  linearly 
interpolated  and  truncated  so  that  is  is  marginally  less  than  the 
time  of  the  concentration.  There  is  also  a  check  to  insure  that  the 
waterway  and  spill  configurations  remain  within  the  limitations  of 
the  microcomputer. 

The  generated  output  and  labelling  information  is  displayed  in  a 
table  which  may  be  reproduced  as  hardcopy  if  desired. 
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Figure  2 


5 . 0    COMFUTAT I DMAL    PROCEDURE 
5.  1    FlQ^^^    Chart    For    Spill     liodel 

M  AME    OF    GUBSTAMœI 


DECAY  RATE 
(Y/M) 


m± 


NAME  OF  WATERWAY 
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I 


STARTING  TIME  ?<  DATE 
EMDINIS  TIME  ?.'.    DATE 

1" 


I  STEADY  STATE  SPILL  (Y/M)!- 


N 


TIME  OF  RELEASE 
MASS  FLUX 


DISCRETE 
COhlCEIMTRATinNS 


It  ÛF  C0NCE^4TRATIaNS 
REQU [ RED 


TIME  TÛ  CDNC. 
DISTANCE  TO  CONG. 


CHECK 
HAS  SPILL  ENDED 


N 


SPILL  TRUNCATION 
SUBROUTINE 


CONCENTRATION 
CALCULATION 
SUBROUTINE 


DECAY  RATE 
SUBROUTINE 
(i  f  appl i  cahl g) 


I 


INTER  DECAY] 
RATE  1 


-3^ 
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ERROR    FUNCTION 
SUBROUTINE 


N 
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I 


PRINT    TABLE     (Y/hÏÏT- 

T 


N 


|PR  I IMT    TAi;<LE    SUBfvOUT  [  ME 


END 
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TIME  OF  RELEASE 
MASS  FLUX 


CONCENTRATION 
PROFILE 


I 


DIST.  TO  PROFILE 


1 


CHECKS; 

IS  SPILL  WITHIN__ 

THE  COMPUTER'S 

CAPACITY 

IS  CONC.  GREATER 

THAN     lE-6 


OK 


Finished 


Aborts 


run 


GOTO; 

OPTION    MENU- 
END— J 
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5.2  Limitations  of  the  Program 

Part  of  the  spill  concentration  is  calculated  by  combining  three 
sub-equations. 

CONC  =  {T2*(T3-T4)  (1) 

T2=EXP(Ve1ocity*Distance/Coeff) 

T3=ERF(SQR(B*T)   +  SQR    (A/T)) 

T4=  ERF(SQR(B*(T-TF))   +  SQR(A/(T-TF  )  )  ) 

Velocity  =  Watercourse  velocity 

Distance  =  Distance  to  Point  of  Concentration  Calculation 

Coeff    =  Dispersion  Coefficient 

B       =  Vel**2/(4*Coeff) 

A       =  Distance**2/(4*Coeff) 

T      =  Time  of  Travel  to  Point  of  Concentration  Calculation 

TF      =  Time  of  Release  of  the  Spill 

ERF     =  Error  Function  Subroutine:   ERF(>2.5)  =  1.0 

The  maximum  real  number  allowed  on  a  micro-computer  before  overflow 
occurs  is  1.7E38.  The  "T2"  calculation  will  result  in  overflow  if 
(Velocity*Distance/Coeff)  exceeds  88  .  T2  will  automatically  be 
assigned  a  value  of  1.7E38,  which  is  incorrect. 

It  has  been  observed  in  many  varying  test  simulations  that  as 
velocity  and  distance  increase  to  cause  overflow,  "T3"  and  "T4" 
become  equal.  Referring  to  equation  1,  if  "T3"  and  "T4"  are  equal 
then  the  value  of  T2  is  irrelevant  as  T2*  (T3-T4)  =  0.  There  is  an 
error  checking  facility  in  the  program  to  insure  if  there  is 
overflow  that  "T3"  and  "T4  are   equal.  This  subroutine  will  not 
insure  that  the  waterway's  configuration  will  always  be  within  the 
computers  limitations,  but  will  increase  the  versatility  of  the 
algorithm  to  include  many  configurations  previously  beyond  the 
capacity  of  the  system. 
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5.3  Assumptions  and  Approximations  Inherent  in  the  Computer  Program 

If  the  argument  of  an  "Error  Function"  exceeds  2.5,  the  program 
assigns  a  value  of  1.0  to  the  returning  error  function  argument. 

If  the  spill  is  not  completed  at  the  time  a  concentration  profile  is 
calculated,  the  remaining  flux  and  time  of  the  spill  is  truncated. 
In  order  to  avoid  division  by  zero,  the  time  of  travel  to  the 
concentration  point  is  increased  in  magnitude  by  l.OE-9  seconds. 
This  addition  to  the  time  of  travel  has  no  effect  upon  the 
concentration  calculation  as  the  calculation  is  performed  at  the  new 
time. 
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6.0  UNITS  OF  MEASUREMENT 

This  program  uses  the  S.I.  system  of  measurement  for  all  input  and 
output  data. 

INPUT  REQUIREMENTS 

Average  Channel  Width   ''  METERS 

Average  Channel  Depth  METERS 

Stream  Velocity  METERS/SECOND 

Dispersion  Coefficient  METERS  SQUARED/SECOND 

Decay  Rate  /SECOND 

Distance  From  Spill  Source  METERS 

Time  Since  Spill  HOURS 

Time  of  Release  of  Spill  HOURS 

Mass  Flux  of  Spill  GRAMS/SECOND 


To  assist  in  the  conversion  of  Imperial  Units  to  the  SI  system  refer 
to  the  following  table. 


MULTIPLY 
BY 

0.0254 

0.305 

1610.0 

0.305 

0.4472 

0.09303 

720.03 

453 

28.313 


FROM 

inches 

metres 

feet 

metres 

miles 

metres 

feet/sec 

m/s 

mph 

m/s 

ft**2/sec 

m**2/s 

mile**2/hr 

m**2/s 

Ibs/s 

g/s 

oz./s 

g/s 
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7.0  ERRORS  AND  ACCURACY 

7.1  Checking  Features  in  the  Program 

All  of  the  input  data  is  checked  to  ensure  that  only  logical  values  can 
be  entered.  The  decay  rate  (if  applicable)  must  be  a  real  number 
greater  or  equal  to  zero.  The  waterway  width,  depth,  velocity  and 
dispersion  coefficient  must  all  be  positive  real  numbers  greater  than 
zero,  if  the  value  entered  does  not  meet  these  criteria  the  question  is 
reprinted. 

The  number  of  mass  flux  changes  entered  must  be  a  positive  integer 
greater  than  zero  and  the  time  of  release  for  each  mass  flux  must  be  a 
positive  real  number  greater  than  zero.  The  mass  flux  must  be  a  real 
number  greater  than  or  equal  to  zero.   If  an  invalid  value  or 
non-numeric  character  is  entered  for  any  of  these  parameters  the 
question  will  automatically  be  reprinted. 

The  distance  downstream,  the  time  of  the  concentration  calculation  and 
the  time  increment  for  the  concentration  profile  must  be  positive  real 
numbers  greater  than  zero  or  the  question  will  be  reprinted. 

The  program  compares  the  spill  time  to  the  time  of  the  concentration 
calculation  to  insure  that  the  spill  has  ended.  If  the  spill  is  still 
in  progress  at  this  point  in  time  the  remainder  of  the  mass  flux  and 
spill  time  is  truncated.   If  the  time  to  the  concentration  calculation 
is  less  than  zero,  the  time  is  changed  to  zero  and  the  profile  is 
considered  completed. 
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7.2  Error  Messages 

If  a  "BEEP"  sounds  and  the  input  parameter  question  reappears  on  the 
screen,  an  illegal  value  or  character  has  been  entered  for  that 
parameter.  Refer  to  section  7.1  of  the  Engineering  and  Software 
Documentation  and  section  3.0  of  the  Users  Manual  for  proper  input  data. 

"SPILL  CLOUD  UNDETECTABLE  (check  input  data)"  If  this  message  appears 
then  the  maximum  concentration  (which  is  assumed  to  occur  at  the  fluid 
travel  time  (distance/velocity)after  the  start  of  the  spill)  does  not 
exceed  lE-6  mg/1 . 

"THIS  CONFIGURATION  EXCEEDS  THE  CAPACITY  OF  THIS  COMPUTER  (refer  to 
users  manual)".  Refer  to  section  5.2  of  this  Manual  for  the  explanation 
as  to  the  cause  of  this  message.  Corrective  measures  are   limited  to 
changing  input  parameters.  The  velocity  and/or  distance  downstream  to 
the  concentration  calculation  can  be  decreased  and/or  the  dispersion 
coefficient  can  be  increased. 

7.3  Numerical  Error 

All  real  number  input  parameters  and  real  number  calculations  are   double 
precision  formatted.  This  causes  all  internal  calculations  to  be 
maintained  at  17  decimal  digits.  The  output  concentration  values  are 
only  displayed  to  three  decimal  digits  (example:  .767D-6).  This  causes 
rounding  error  and  possible  variation  between  identical  simulations. 

Due  to  the  approximations  involved  with  integrating  the  Fickian  spill 
equation  to  account  for  a  finite  time  release  spill,  results  obtained 
from  identical  simulations  may  not  exactly  match.  This  affect  may  be 
amplified  when  identical  simulations  are   calculated  on  different  models 
or  brands  of  microcomputers  depending  upon  the  logic  and  accuracy 
programmed  into  the  computer. 
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7.4  Suggestions  to  Improve  the  Accuracy  of  the  Simulation 

Since  the  model  relies  entirely  on  input  data  to  operate,  input 
parameters  that  are   accurate  and  precise  can  greatly  improve  the 
quality  of  the  simulation.   It  is  recommended  that  all  data  be 
entered  with  a  sufficient  number  of  significant  digits  to  reflect 
the  double  precision  format  of  the  internal  calculations  (example: 
desirable:  1.3276E3,  undesirable:  1.3E3). 

It  is  recommended  that  the  spill  be  divided  into  small  enough 
segments  such  that  every  segment  is  essentially  steady  state. 

Physical  stream  parameters  (width,  depth,  velocity  and  dispersion 
coefficient)  should  reflect  the  average  characteristics  of  the 
entire  reach  from  spill  source  to  the  point  of  the  furthest 
concentration  calculation  downstream. 

A  dispersion  coefficient  must  be  derived  analytically  if  there  is  no 
experimental  simulation  data  available  for  the  reach.  A  suggested 

5 

equation  to  calculate  a  dispersion  coefficient  is: 

COEFF  =  5*U*W**2/d 

Coeff:  Dispersion  Coefficient 

U:  Shear  Velocity  (approx.=  0.1*  Stream  Velocity) 

W:  reach  width 

d:  reach  depth 

If  there  is  no  background  information  to  base  an  estimate  for  a 
dispersion  coefficient  on,  an  initial  value  of  5.0  m**2/s  should  be 
used.  Turbulent,  fast  flowing  rivers  have  larger  coefficients 
(10-30  m**2/s).  Slow  moving,  relatively  straight  reaches  usually 
have  lower  coefficients  (1-5  m**2/s). 

If  a  decay  rate  is  utilized,  a  safety  factor  should  be  introduced  by 
underestimating  the  decay  rate.  The  minimum  recommended  decay  rate 
should  always  be  utilized. 
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7.5  Accuracy  of  the  Model 

This  model  was  executed  with  the  same  input  data  as  two  simulations 
(Ottawa  and  Speed  Rivers)  that  were  run  with  the  FORTRAN  model  on 
the  mainframe  computer.  The  micro-computer  version  duplicated  the 
mainframe  results  exactly  (figures  3,4)  .  The  original  mainframe 
results  for  the  Speed  River  was  verified  against  observed  data 
(figure  5). 
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Figure  3 

Comparison  of   the  Mainframe  Versus  Micromputer  Results  for  the 

Speed  River 
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Figure  4 

Comparison  of   the  Mainframe  Versus  Microcomputer  Results  for  the 
Ottawa  River 
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Figure  5 

Concentration  Profile  for  Simulated  and  Observed  Dye  Test  Data  in  the 

Speed  River' 
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8.0  SOFTWARE  AND  HARDWARE  REQUIREMENTS 

The  program  is  supplied  on  a  5  1/4  inch  IBM  formatted  floppy  disk. 
The  compiled,  ready  to  run,  version  of  the  program  and  a  program 
1 isting  are  suppl ied. 

Editing  the  program  listing  will  have  no  effect  upon  the  compiled 
version.  The  altered  program  listing  must  be  recompiled  in  order 
for  changes  to  be  effective.  The  accuracy  and  reliability  of  this 
product  cannot  be  guaranteed  if  it  is  altered  or  changed.  The  user 
should  recheck  the  program  with  the  test  data  supplied  in  the  Users 
Manual  after  editing  the  model. 

The  hardware  required  for  this  model  is  an  IBM  or  IBM  compatible 
micro-computer,  running  MS  DOS,  with  at  least  256  K  of  memory.  A 
printer  is  necessary  if  a  hardcopy  of  the  displayed  output  is 
required. 
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APPENDIX  1.0 

Lateral  Mixing 

The  longitudinal  distance  to  the  point  of  complete  natural  lateral 
mixing  depends  upon  such  factors  as  reach  velocity,  depth,  width, 
turbulence,  and  relative  straightness  of  the  channel.  A  suggested 
equation  to  calculate  the  distance  to  this  point  is: 

Lateral  Mixing  Distance  =  (V*W**2)/Lat  Coeff 

V  =  Average  Reach  Velocity 

W  =  Average  Reach  Velocity 

Lat  Coeff  =  A*d*(Shear  V) 

d  =  Average  Reach  Depth 

Shear  V  =  (0.1*Average  Reach  Velocity) 

A  =  .1  to  .2  for  Straight  Channels 

.6  For  Most  Streams 

(irregularities  cause  the  coefficient  to  increase) 


Example 

Velocity  =  .1  m/s 

Width  =  30  m 

Shear  Velocity  =  (.l*.l  m/s)  =  .01  m/s 

Depth  =  2.0  m 

A  =  .6 

Lat  Mixing  Dist  =  ( .l*30**2)/( .6*2.0*. 01 ) 
=7500  m 

In  this  example,  the  one  dimensional  spill  model  can  be  considered 
valid  7500  m  downstream  from  the  spill  source. 
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Program  Listing 
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APPENDIX     2.0 

Program  Listing 
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10  PBINT  ■lllllllllllllllllllllllIllllllllllllllllllllllllllllllllHllllllllllll 

111111111" 

20  PBINT  ■llllllllllllllllllllllllllllllUllllllllllllllllllllllllllllllllllllll 

lllllllir 

30  PRINT  "1SSSSSS3SSSSSSSSSSS1PPPPPPPPPPPPPPP1IIIIII1LLLLLLLLLLLLLLL1LLLLLLLLLLL 

LLLLLLLLl" 

PlI   IlL   LLLLLLLLLLIL   LLLLLL 


SIP 
SIP 


PlI   IlL   LLLLLLLLLLIL   LLLLLL 


SIP 


PlI   IlL   LLLLLLLLLLIL   LLLLLL 
PlI   IIL   LLLLLLLLLLIL   LLLLLL 


SIP   PPPPPPPPPPII   IlL   LLLLLLLLLLIL   LLLLL 


40  PRINT  'ISSSSSS 

LLLLLLLLl" 

50  PBINT  "ISSSSSS 

LLLLLLLLl' 

60  PBINT  "ISSSSSS   SSSSSSSSIP   PPPPPP   PlI   IlL   LLLLLLLLLLIL   LLLLLL 

LLLLLLLLl" 

70  PRINT  "ISSSSSS   SSSSSSSSIP   PPPPPP   PlI   IlL   LLLLLLLLLLIL   LLLLLL 

LLLLLLLLl" 

80  PBINT  "ISSSSSS   SSSSSSSSIP 

LLLLLLLLl" 

90  PRINT  "ISSSSSS 

LLLLLLLLl" 

100  PBINT  "ISSSSSS 

LLLLLLLLLl" 

110  PRINT  "ISSSSSSSSSSSSS   SIP   PPPPPPPPPPII   IlL   LLLLLLLLLLIL   LLLLL 

LLLLLLLLLl" 

120  PRINT  'ISSSSSSSSSSSSS   SIP   PPPPPPPPPPII   IlL   LLLLLLLLLLIL   LLLLL 

LLLLLLLLLl" 

130  PRINT"1SSSSSSSSSSSSS   SIP   PPPPPPPPPPII   IlL   LLLLLLLLLLIL   LLLLLL 

LLLLLLLLl' 

140  PRINT'ISSSSSS 

LLLLLl" 
150  PBINT'ISSSSSS 

LLLLLl" 
160  PBINT"1SSSSSSSSSSSSSS3SSS1PPPPPPPPPPPPPPP1IIIIII1LLLLLLLLLLLLLLL1LLLLLLLLLLL 
LLLLLLLLl" 

170  PBINT'llllllllllllllllllllllllllllllllllllllllUllllllllllllllllllUllllllll 
111111111" 
180  PBINT 

190  PBINT  "Developed  for  P.C.  by  Ï.C.   AlliBon,  Waterloo  Co-op  Student,  River  Syg 
tens  Unit" 

200  PRINT  "Ontario  Ministry  of  the  Environnent.  April  1986  ' 
210  PBINT  "Ail  enquiries  (in  writing)  should  be  directed  to  Z.  NOVAR" 
220  PRINT  "  Biver  Systens  Unit,  Water  Besources  Branch  " 

230  PBINT  "  Ministry  of  the  Environnent" 

240  PBINT  "  135  St.  Clair  Avenue  West,  Toronto,  Ontario" 

250  INPUT  "HIT  BBTUBN",ASD$ 
260  'tttntttnttsPILL  HODBLnntttnutttt 
270  'ONB  DIMENSION  FINITE  TIME  BBLBASB 

280  '  BASED  ON  THEOBY  OF  SUPEBPOSITION  AND  THE  FICRIAN  DISPEBSION  THBOBY 
290  '  DEVELOPED  JANUARY  to  HAY  1986 
300  '    BY  W.C.  ALLISON,  WATBBLOO  CO-OP  STUDENT 

310  PBINT  "tttmtnuutttJWELCOHB  TO  THE  WORLD  OF  SPILL  HODBLLING«»»t»t»nn 
t" 


SIP   PPPPPPPPPPII    IlL 
SIP   PPPPPPPPPPII    IlL 


LIL 
LIL 
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320  PBINT:PRINT:PSIIIT 

330  PRINT  "This  is  a  one  diiensional  finite  tiie  release  spill  «odel  based  on  th 

e  fickian  dispersion  equations  and  the  theory  of  superposition.' 

340  'define  types  of  variables  and  use  double  precision 

350  DBFDBL  A-Z 

360  DBPINT  P,P,N,I,J 

370  'define  fluï  eiîe  arrays 

380  DIM  THASS(100),TF(100),TFTEHP(100),THTBHP(100) 

390  PKINT:P8INT:PRINT:PBINT:PRINT:PBmT:P8INT:P8INT:PRINT 

400  PBINT:P8mT:PEINT: PRINT: PRINT: PRINT: PRINT 

410  LINE  INPUT  'BNTBB  THE  NAHB  OF  THE  SUBSTANCE  THAT  HAS  BEEN  SPILLED 

(typical  exaiple:  perchloroethelyene)  ';  CHBHJ 
420  '  after  correction  go  back  to  table 
430  IF  NLUl  THEN  GOTO  660 
440  PBINT 

450  'flag  for  decay  rate  application 
460  DFLAG:0 

470  INPUT  'IF  THIS  IS  A  NON-CONSBBVATIVE  SUBSTANCE  DO  YOU  KISH  TO  APPLÏ  AN  BÏPON 
BNTIAL    DECAY  BATB,  if  conservative  type  'N'  (Y/N)",  DEC} 
480  PBINT 

490  IF  (DBC$:'Y*  OB  DBC):"y")  THEN  DFLAGrl  :GOTO  520 
500  IF  (DEO^'N"  OB  DBCJ:"n")  THEN  GOTO  560 
510  BEBP:GOTO  470 

520  INPUT  '«HAT  IS  THE  DECAY  BATB  (s'-l),  SEE  USEES  MANUAL  FOB  BBFBEBNCE';DBCAY 
530  IP  DBCAT<0  THEN  BEEP  :GOTO  510 
540  PBINT 

550  '  after  correction  go  back  to  table 
560  IF  NL5:1  THEN  GOTO  660 
570  LINK  INPUT  "BNTEB  THE  NAME  OF  THE  BODY  OF  WATEB  THB  SPILL  HAS  OCCUBED  IN 

(typical  eranple:  Huaber  River]  ";  RIVEEJ 
580  '  after  correction  go  back  to  table 
590  IF  NL2:1  THEN  GOTO  660 
600  PRINT 
610  LINE  INPUT  'ENTER  THB  NAME  OF  THB  SOURCE  OF  THB  SPILL 

(typical  exaaple:  Rolphton  Nuclear  Plant)  ';  SOURCES 
620  '  after  correction  go  back  to  table 
630  IF  (JL3:1  THEN  PBINT:GOTO  660 
640  PBINT:PBINT:PRINT:PBINT 
650  'table  of  preliminary  data 
660  PBINT  'IITHE  SUBSTANCE  SPILLED  «AS  "CHBH$ 
670  PRINT 

680  PRINT  ■2)THB  NAHB  OF  THB  RIVER  IS  "RIVERJ 
690  PRINT 

700  PRINT  ■3)THB  SPILL  OCCURBD  AT  "SOUBCE$ 
710  PBINT 

720  'bypass  if  no  decay  rate 
730  IF  DFLAG:0  THEN  GOTO  780 

740  PBINT  "4)THB  BÏPONENTIAL  DECAY  BATB  IS  "DBCAY  'S'-l* 
750  PBINT 

760  INPUT  'DO  YOU  «ISH  TO  CHANGE  1,2,3,4  or  5(G0  0N|  ';  Hi 
770  GOTO  790 

780  INPUT  "DO  YOU  WISH  TO  CHANGE  1,2,3,  or  4(G0  ON)  ";  Z2$ 
790  PRINT:PRINT 
800  'directions  to  corrections 


810  IF  (Z2J:"1')  THEN  NLUl 

820  IF  (Z2$:'2")  THEN  NL2:1 

830  IF  (Z2|:'3')  THEN  NL3:1 

840  'decay 


GOTO  410 
GOTO  570 
GOTO  610 
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850  IF  {Ui-.'r   AND  DFLAG:1)  THEN  NL5:1:  GOTO  460 

860  'no  decay 

870  IF  [Ui-.'r   AND  DFLAGOl)  THBN  GOTO  930 

880  'decay 

890  IP  (Z2$:"5"  AND  DFLAG:I)  THBN  GOTO  930 

900  BEEP 

910  IF  DFLAG:!  THEN  GOTO  760 

920  IF  DFLAGiO  THEN  GOTO  780 

930  PRINT:P8INT:PBINT:PBINT 

940  PRINT  "THE  SPILL  START  AND  END  TIMES  AND  DATES  HUST  BE  ENTERED" 

950  PSINT:PRINT:PRINT 

960  LINE  INPUT  "ENTER  STARTING  TIHE  AND  DATE  OF  SPILL  (typical  exaiple;  2:30  pa 

January  11,1985)";  STARTJ  :IF  (FLUl)  THBN  GOTO  990 

970  LINE  INPUT  "ENTER  ENDING  TIHE  AND  DATE  OF  SPILL  (typical  exanple;  8:30  pa  Ja 

nuary  11,1985)  ";  SPEND$ 

980  PRINT:PRINT:PRINT 

990  PRINT  "THE  SPILL  LASTED  FROM" 

1000  PRINT  TAB(20)  START} 

1010  PRINT  "  TO' 

1020  PRINT  TAB(20)  SPEND) 

1030  INPUT  'DO  YOU  KISH  TO  CHANGE  EITHER  OF  THESE  TIHBS  AND  DATES  (Y/N)";  Zl$ 

1040  PRINT:PRINT:PRINT 

1050  IF  (Z1J:"Y"  OR  Zl$:"y*)  THEN  GOTO  1080 

1060  IP  (Z1$:"N"  OR  m-.'n']   THEN  GOTO  1140 

1070  BEEP:  GOTO  1030 

1080  INPUT  'DO  YOU  WISH  TO  CHANGE  1)  START  2)BND  3)N0  CHANCE";  Z$ 

1090  PRINT:PRINT:PRINT:PEINT:PRINT:PRINT 

1100  IF  [li-'V]   THEN  FLUl:  GOTO  960 

1110  IF  (Z$:"2")  THEN  PL2:1:  GOTO  970 

1120  IF  (Z$:"3")  THEN  GOTO  1140 

1130  BBBP:GOTO  1080 

1140  PRINT:PRINT:PRINT:PRINT: PRINT: PRINT 

1150  'flags  for  streaa  paraaeter  corrections 

1160  PL1:0:PL2:0:PL3:0:PL4:0 

1170  PRINT  "YOU  HUST  NOW  ENTER  THE  STEBAH  CHARACTERISTICS" 

1180  PRINT:PRINT:PRINT 

1190  INPUT  "THE  CHANNEL  WIDTH  (in  aeters)  IS";  W 

1200  IF  (W<:0)  THEN  BEEP  ;GOTO  1190 

1210  IF  (PLUlj  GOTO  1350 

1220  PRINT:PRINT:PR1NT 

1230  INPUT  "THE  CHANNEL  DEPTH  (in  aeters)  IS";  D 

1240  IF  (D<:0)  THBN  BBEP  :GOTO  1230 

1250  IF  (PL2rl)  GOTO  1350 

1260  PRINT:PRINT:PRINT 

1270  INPUT  "THE  VELOCITY  (in  «eters/second)";  VBL 

1280  PRINT:PRINT:PE1NT 

1290  IP  (VBL<:0)  THBN  BEEP  :GOT0  1270 

1300  IF  (PL3:1)  GOTO  1350 

1310  INPUT  "THE  DISPERSION  COEFFICIENT  (in  aeters  squared  per  second)  IS";  COBFF 
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1320  PEDiT:PEINT:PHINT 

1330  IF  (COEFF<:0|  THEN  3BEP  :GOTO  1310 

1340  IF  (PL4:1)  GOTO  1350 

1350  PBINTiPEINT 

1360  'table  for  streaa  paraaeters 

13?0  PRINT  'STBEAH  CHAEACTERISTICSttttJm" 

1380  PRINT  "DTHB  CHANNEL  IS"  Ï  "aeters  ÏIDB" 

1390  PRINT  "2)THE  CBANNBL  IS"  D  "aeters  DEEP" 

1400  PRINT  "31THE  VELOCITY  OF  THE  STBEAH  IS"  VBL  'leters  per  second" 

1410  PRINT  "4)THE  DISPERSION  COEFFICIENT  OF  THE  STREAM  IS"  COSFF  "aeters  squared 

per  second" 

1420  PBINT:PEINT:PBINT 

1430  INPUT  "DO  TOU  ÏISH  TO  CHANCE  AN7  OF  THESE  PARAHETBRS  (Y/N)";  PU 

1440  PRIMT:PRmT:PRINT 

1450  IP  (P1$:"Y"  OB  PU:"y")  THEN  GOTO  1480 

1460  IP  (P1$:"H'  OB  Pl$:"n")  THEN  GOTO  15?0 

1470  BBBP:  GOTO  1430 

1480  INPUT  "XBICH  PABAHETBR  DO  YOU  ÏISH  TO  CHANGE  (1,2,3,4  or  5(No  Change) )";P2J 

1490  PRINT:PEINT:PRINT 
1500  IP  (P2$:"r)  THEN  PLUl:  GOTO  1190 
1510  IP  (P2$:"2")  THEN  PL2:1:  GOTO  1230 
1520  IP  (P2$:"3")  TBBH  PL3il:  GOTO  1270 
1530  IP  (P2$:"4")  THEN  PL4:1:  GOTO  1310 
1540  IF  (P2$:"5")  GOTO  1570 
1550  BEEP:GOTO  1480 

1560  'CALCULATE  THE  STEBAH  ABBA  AFTER  STBEAH  CHARACTBEISTICS  ABB  FINALIZED 
1570  ARBA:li»D 

1580  PRINT  "THE  SPILL  CHARACTERISTICS  HUST  NOK  BE  ENTERED" 
1590  PBINT:PRINT:PRINT 
1600  '  shortcut  if  steady  state 

1610  INPUT  "VAS  THIS  A  STEADY  STATE  SPILL  (ie.  was  there  no  variation  in  the  aas 
8  flux  of   the  substance  spilled)  (Y/HI";  P6$ 
1620  '  steady  state  flag 
1630  SPPLAG:0 

1640  IP  (P6$:'Y"  OR  P6|:"y")  THEN  N:1:SPFLAG:1  :  GOTO  1700 
1650  IP  (P6J:"N"  OB  P6$:"n")  THEN  GOTO  1670 
1660  BEEP  :GOTO  1610 
1670  PRINT:PBINT:PRINT 
1680  '  non-linear  spill  input 

1690  INPUT  "HOU  MANY  MASS  FLUX  CHANGES  WHERE  THERE  (including  interaediate  fluïe 
8  of  zero   graas  per  second)  DURING  THE  SPILL";  N 
1700  PBINT:PEINT;PBINT 
1710  IF  (N<:0)  THEN  BBBP  :GOTO  1690 

1720  INPUT  "BEHBHBBB  ALL  MASS  FLUXES  HUST  BE  SPECIFIED  IN  (graas  per  second)  AND 
HASS  FLUI  DUEATIONS  HUST  SPECIFIED  IN  (hours)  (eiaaple;  56.4  g/s  and  2.85  hour 
s)       HIT  BBTUBN  tumnimt'',P3$ 
1730  'entering  spills 
1740  '  counter  for  spill  input 
1750  FOR  1:1  TO  N 
1760  PRINT:PRINT 

1770  PRINT  'MUtnitFOR  HASS  SPILL  NUHBBRtJnttnn",  I 
1780  PB1NT:PRINT 

1790  INPUT  "THE  DURATION  OP  THIS  MASS  SPILL  IN  HOURS  WAS   ",  TP(I) 
1800  PRINT:PRINT 

1810  IF  (TF(I)(:0)  THEN  BBBP  :PRINT  "REENTER  HASS  SPILL  TIHB  FOR  NUHBBSM  : 
GOTO  1790 
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1820  PRINTiPBIMT 

1830  INPUT  'THB  MASS  FLUÏ  OF  THIS  PORTION  OF  THE  SPILL  (in  grass  per  aecondl  «AS 

",THA3S(I| 
1840  PRINT:P8INT 

1850  IF  (THASS(I)<0)  THEN  BBBP:  PRINT  "RBBNTBR  HASS  FLUX  FOR  NUMBER",!  :GOTO  183 
0 

I860  PBINT:PRINT 
1870  NBXT  I 
1880  'spill  table 

1890  PRINT  "SPILL       TIME  LENGTH       HASS  FLU!" 
1900  PRINT  "NUMBER        hours  ils      " 

1910  PRINT 
1920  FOB  1:1  TO  N 

1930  PRINT  I  ,'  "TFdl,"      'THASS(I) 
1940  NBÏT  I 

1950  INPUT  "DO  YOU  WISH  TO  ADD  TO  OR  CHANGE  ANY  OF  THESE  VALUES  (Y/N)   ';  ïl$ 
1960  P8INT:PRINT:PRINT 
1970  IF  (Ï1$:"Y"  OR  ïlj:"y")  THEN  GOTO  2010 
1980  IF  (ïl$:'N"  OB  n$:"n"i  THEN  GOTO  2190 
1990  BEEP  :GOTO  1950 
2000  PRINT:PRINT 

2010  INPUT  "YOU  MAY  1)ADD  ANOTHER  SPILL  or  2)CHANGB  ONE  or  3|G0  ON    ',  Ui 
2020  PRINT:PRINT:PRINT 

2030  IF  (X2$:°r)  THEN  PRINT  "YOU  MUST  REENTER  ALL  SPILLS" :GOTO  1690 
2040  IP  (Ï2$:"r  AND  SPFLAG:1)  THEN  11=1  :GOTO  2110 
2050  IF  (Ï2$:"2*)  THEN  GOTO  2080 
2060  IF  (Ï2}:"3')  THEN  GOTO  2190 
2070  BBBPiGOTO  2010 

2080  INPUT  "ÏHAT  IS  THB  NUMBER  OF  THE  SPILL  YOU  )»ISH  TO  CHANGE  ";I1 
2090  PRINTiPRINT 

2100  IF  (I1(:0)  OR  (I1>N)  THEN  BEEP  :GOTO  2080 
2110  INPUT  "THE  NEW  DURATION  OF  THIS  SPILL  IS  (in  hours)  ',  TF(Il) 
2120  PEINT:PRINT 

2130  IP  (TF(I1)<:0)  THEN  BBBP  :GOTO  2110 
2140  INPUT  "THE  NEW  HASS  FLUX  OP  THIS  PORTION  OF  THB  SPILL  (in  grans  per  second) 

IS  ",  TMASSdD 
2150  PRINT:PR[NT:PRINT 
2160  IF  (THASS(I1|<0)  THEN  BBBP  :GOTO  2140 
2170  GOTO  1890 

2180  'superiupose  tine  of  spill  fluxes  and  convert  to  seconds 
2190  FOR  1:1  TO  N 
2200  TF(0):0 

2210  TP(I):TF(I)«3600tTF(I-l) 
2215  TPTEMP(I):TP(I)  :  THTBMP(I):THASS(I) 
2220  NBXT  I 

2230  'DIMENSION  REMAINING  ARRAYS 

2240  DIM  T(200),DIST(200),TCONC(200),TCON(2Û0),Tl(200),T2(100),T3(100),TCONCl(20 
0|,TCONC2(IOO),TCONC3(100),TT(200),TF1(100)JMASS1(100) 

2250  'GENERAL  RULES  FOR  VARIABLES;  if  there  is  a  "c"  in  it  concentration,  "tf"-t 
iae  of  spill  fluxes,  '«"-spill  lass  fluxes,  "t'-alone-tine  to  concentration,  "te 
mp-teiporary  storage  array 
2260  INPUT  "THB  REQUIRED  INPUT  DATA  IS  NOW  COMPLETE,  YOU  MUST  NOW  DECIDE  ON  WHAT 

TYPE  OF  OUTPUT  YOU  REQUIRE  FROM  THB  OPTIONS  THAT  WILL  APPEAR  ON  THB  SCREEN. 

HIT  RETURN!U«tJnimunut  " ,m 
2270  PRINTiPRINT 

2280  PRINT  "tttmOPTIONStttm" 
2290  PRINTiPRINT 

2300  PRINT  "IIDISCRETE  CONCENTRATIONS  FOR  A  SERIES  OF  DISTANCES  ANU  TIMES  DOWNST 
BEAM" 
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2310  PRINT  "2|A  CONCENTBATION  PROFILE  OVER  TIHE  FOR  A  DISTANCE  DOWNSTREAM" 

2320  PRINT 

2330  INPUT  "DO  YOU  KISH  OPTION  1  OR  2";  VU 

2340  IF  (Vl}:"!")  THEN  GOTO  2370 

2350  IP  (Vl$:"2")  THEN  GOTO  3020 

2360  BEEP  :GOTO  2280 

2370  PRINT:PRINT 

2380  'discrete  concentrations 

2390  PRINT  "ALL  TINES  HUST  BE  SPECIFIED  RELATIVE  TO  THE  STARTING  TIHE  OF  THE  SPI 

LL.  (exaaple;  specifying  a  tine  of  60  hours  and  5000  neters  would  produce  a  cone 

entration   5000  neters  downstreaa  froii  the  spill,  60  hours  after  the  spill  star 

ted" 

2400  PRINTlPRINT 

2410  INPUT  "HOW  MANY  DOWNSTREAM  CONCENTRATIONS  DO  YOU  DESIRE  ";  NTC 

2420  PRINTlPRINT 

2430  IF  (NTC(:0)  THEN  BEEP  :GOTO  2390 

2440  'counter  for  nuaber  of  concentrations 

2450  FOR  Irl  TO  NTC 

2460  PRINT  "FOR  CONCENTRATION  NUHBERnitniinu  ",[ 

2470  PRI(IT:PRINT 

2480  INPUT  "THE  TIHE  (in  hours)  SINCE  THE  SPILL  STARTED  IS  ";T(I1 

2490  PRINT:PRINT 

2500  IF  (T(I)<:0)  THEN  BEEP  :GOTO  2460 

2510  INPUT  "THE  DISTANCE  DOWNSTREAM  TO  THE  PLACE  OF  THE  DESIRED  CONCENTRATION  (i 

n  seters)  IS   ";DIST(I) 

2520  IP  (DIST(I)<:0)  THEN  BEEP  :GOTO  2510 

2530  NEIT  I 

2540  PRINT:PRINT:PRINT 

2550  'table  for  desired  discrete  concentrations 

2560  PRINT  'NUMBER      TIME       DISTANCE 

2570  PRINT 

2580  FOR  1:1  TO  NTC 

2590  PRINT  I,  "   "T(I),  "    "DISTd) 

2600  NEXT  I 

2610  INPUT  "DO  YOU  WISH  TO  1)CHANGE  A  VALUE  21ADD  ANOTHER  TIHE  3|G0  ON  ";V21 

2620  PRINT 

2630  IF  (V2$:"l')  THEN  GOTO  2630 

2640  IF  (V2$:"2")  THEN  PRINT  "THE  ENTIRE  LIST  MUST  BE  RETYPED"  :GOTO  2390 

2650  IP  (V2$:'3")  THEN  GOTO  2800 

2660  BEBPiGOTO  2610 

2670  PRINTlPRINT 

2680  INPUT  "WHAT  TIME  NUMBER  DO  YOU  WISH  TO  CHANGE  ";I1 

2690  PRINTlPRINT 

2700  IF  (I1<:0)  OR  (IDNTCI  THEN  BEEP  :GOTO  2680 

2710  PRINT  "FOR  CONCENTRATION  NUMBERntJtnnmn"  H 

2720  PRINTlPRINT 

2730  INPUT  "THE  NEW  TIME  TO  THE  REQUIRED  CONCENTRATION  IS   ",  T(I11 

2740  PRINTlPRINT 

2750  IF  (T(I1)<:0)  THEN  BEEP  iGOTO  2710 

2760  INPUT  "THE  NEW  DISTANCE  TO  THE  REQUIRED  CONCENTRATION  IS   ",DIST(11)  , 

2770  PRINTlPRINT 

2780  IF  (DIST(I11<:0)  THEN  BEEP  iGOTO  2760 

2790  GOTO  2560 
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2800  PBINTlPRINT 

2810  'go  to  concentration  subroutine 

2820  COSUB  3760 

2830  'convert  tiae  and  concentrations  into  final  arrays  to  be  used  in  output 

2840  FOB  1:1  TO  NTC 

2850  TCON(I):TCONC(I) 

2860  TT(I):T(I)/3600 

2870  NEXT  I 

2880  NJrNTC 

2890  'flag  for  distance  to  be  in  table  for  discrete  output 

2900  FLAGrl 

2910  'go  to  output  table  subroutine 

2920  GOSUB  4700 

2930  •  reset  flag  to  0 

2940  FLAG:0 

2950  PEINT:PEINT 

2960  INPUT  "DO  YOU  UANT  1)M0RB  DISCRBTB  CONCENTRATIONS  2)EBTURN  TO  OPTION  HENU  3 

)QUIT",X3$ 

2970  PRINT:PRINT 

2980  IF  (Ï3<r'r|  THEN  GOTO  2390 

2990  IF  (X3$:'2''|  THEN  GOTO  2260 

3000  IF  (X3$:"3")  THEN  GOTO  3750 

3010  BEEP  :GOTO  2960 

3020  'concentration  profile  option 

3030  PRINT:PEINT:PRINT 

3040  '  position  of  profile 

3050  INPUT  "WHAT  DISTANCE  DOWNSTREAM  (in  leters)  DO  YOU  WISH  TO  PROFILE  OVER  TIH 

E";  DISTdl 

3060  PRINTlPRINT 

3070  IF  (DIST(1)<:0)  THEN  BEEP  :GOTO  3050 

3080  INPUT  "WHAT  TIME  STEP  BETWEEN  CONCENTRATIONS  (in  hours)  DO  YOU  WISH,  Large 

spills  should  have  larger  tine  steps";  TS 

3090  IF  (TS<:0)  THEN  BEEP:  GOTO  3080 

3100  'using  velocity  and  distance  to  profile  find  an  initial  starting  position  d 

ownstrean,  the  prograa  works  in  both  directions  fron  this  point 

3110  NTC:1 

3120  T(1):DIST(1)/VEL/3600 

3130  'SET  INITIAL  VALUES  ASIDE 

3140  TKBEPiT(l) 

3150  GOSUB  3760 

3160  CKBBP:TC0NC(1) 

3170  'SET  «INIHUM  DETECTABLE  CONCENTRATION 

3180  FLAGT:0 

3190  IF  (TCONCdl):. 000001)  THEN  GOTO  3240 

3200  BEEP 

3210  PRINT  "SPILL  CLOUD  UNDETECTABLE  TO  l.OE-6  ig/1  (check  input  data) 

32^20  GOTO  2960 

3230  'increaent  forward  over  tiae  froa  initial  value 

3240  12:1 

3250  T(1):TKBEP:3600 

3260  T(l):{T(l)/360û)tTS 

3270  N2:I2 

3280  'go  to  concentration  subroutine 

3290  GOSUB  3760 

3300  T2(I2):T(1) 
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3310  TCONC2(I2):TCONC(l) 

3320  I2:I2tl 

3330  'repeat  until  cloud  ends 

3340  IP  (TCONCd)):. 000001)  THEN  GOTO  3260 

3350  13:0 

33B0  T(1):TRBEP«3600 

3370  'increnent  back  in  tiae  fron  initial  value 

3380  T(lj:(T(l|/360û)-TS 

3390  'if  tine  becoaes  negative  set  to  jero 

3400  IP  T(1)<0  THEN  T(1):0 

3410  I3=I3tl 

3420  'goto  concentration  subroutine 

3430  GOSUB  3760 

3440  T3(I3):T(1| 

3450  N3:I3 

3460  TCONC3(I3):TCONC(l) 

3470  'repeat  until  cloud  disappears 

3480  IP  (TC0NC(1)>:. 000001)  THEN  GOTO  3380 

3490  'set  up  final  array  for  table  output 

3500  I5:M3 

3510  14:1 

3520  KHILE  (I5>:1)  AND  (I4<:N3) 

3530  TT(I4):T3(I5) 

3540  TT(I4):TT(I4)/3600 

3550  TCON(I4):TCONC3(I5) 

3560  14:14+1 

3570  15:15-1 

3580  GOTO  3520 

3590  «END 

3600  'set  initial  value  in  array 

3610  TT(N3tl):TKBBP  :TC0N(N3U):CRBEP 

3620  FOa  16:1  TO  U 

3630  T2(I6):T2(I6)/3600 

3640  TT(N3  +  UI6|:T2(I6) 

3650  TCON(N3H  +  I6):TCONC2(I6) 

3660  NEXT  16 

3670  NJ:N3+1+N2 

3680  'goto  output  subroutine 

3690  GOSUB  4700 

3700  PRINT:PRINT:PEINT 

3710  INPUT  "DO  YOU  ÏISH  TO  1)G0  TO  OPTION  MENU  2)QUIT";  Z10$ 

3720  IF  (Z10$:"l")  THEN  GOTO  2260 

3730  IF  (Z10J:"2")  THEN  GOTO  3750 

3740  BEEP:GOTO  3710 

3750  END 

3760  'SUBEOUTINB  P08  CALCULATING  CONCENTRATION 

3770  'flag  for  inconplete  spills 

3780  TBUNC:0 

3790  'store  original  nuaber  of  spill  fluxes 

3800  NQ:N 

3810  '  store  original  spill 

3850  PEINT 

3860  PRINT'CALCULATING  CONCENTRATIONS  PLEASE  WAIT" 

3870  PRINT 

3880  'counter  for  nuaber  of  concentrations  required  (1  if  profile) 

3890  FOR  lUl  TO  NTC 

3900  '  convert  to  seconds 

3910  T(I1):T(I1)»3600 

3920  'initialiie  concentration 
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3930  TCONC(I1):0 

33*0  'reassign  original  n 

3950  M-KQ 

3960  'reassign  original  spill  for  each  run  thru 

3970  FOR  112:1  TO  N 

3980  TF(I12):TFTB«P(n2)  :THâSS(I121:THTEHP(112) 

3990  NBÏT  112 

4000  FOB  1:1  TO  N 

4010  'if  spill  is  incoiplete  goto  truncation  subroutine 

4020  IF  T(I1)<:TF(I)  THEN  GOSUB  5470 

4030  'initial  calcs 

<0<0  A^(DrST(Iir21/(4lCOBFF) 

4050  B:(VBL*2)/(4»C0BFF) 

4060  T1:TMASS(I)/(2»ABBA«VBL) 

4070  VDC:VBLJDIST(I1)/C0BFF 

4080  'PC  liiitations  allows  a  sax  of  exp{88)  see  users  lanual  for  bypass  aethod 

4090  IP  VD088  THBN  T2:0  :  GOTO  4110 

4100  T2:BîP(VDC) 

4110  T3:SQR(BtT(Ill)+SQE(A/T(Il|) 

4120  Ï:T3 

4130  'goto  error  function  subroutine 

4140  GOSUB  4330 

4150  T3:BRF 

4160  T4:SQB(BJ(T(I11-TF(I1))+SQR(A/(T(I1)-TF(I1)) 

4170  I:T4 

4180  GOSUB  4380 

4130  T4iBEP 

4200  'overflow  protection  see  users  iianuai 

4210  IF  ((T3-T4)<>0  AND  VDC>88)TBBN  FEINT  "THIS  CONFIGURATION  BXCBBDS  THB  CAPACI 

T!  OF  THIS  COHPUTBR  (refer  to  users  aanuai)'  :  GOTO  3750 

4220  T5:SQE(B»T(I1))-SI3E(A/T(I1!) 

4230  ï:T5 

4240  GOSUB  4380 

4250  T5:BBF 

4260  T6:SQB(BJ(T(I1)-TF(I)))-SQR(A/(T(I1)-TF(I)1) 

4270  ï:T6 

4280  GOSUB  4380 

4290  T6:BRF 

4300  'sui  up  to  get  concentration 

4310  CONC:Tlt(T2J(T3-T41+T5-T6) 

4320  'sus  up  concentrations  fros  each  flux  and  set  in  array 

4330  TC0NC(Il):TCONC(Il)tCONC 

4340  NBÏT  I 

4350  NEXT  II 

4360  RETURN 

4370  ' 

4380  'subroutine  for  error  function 

4390  'source  is  referenced  in  iianual 

4400  l-l 

4410  NBG:0 

4420  IF  (Z>:0)  THBN  GOTO  4460 

4430  Z:Zt(-l) 

4440  NBG:1 

4450  'if  input  f  is  >2.5  erf  goes  to  I  m  iniinity  so  we  issign  it  1  at  this  poi 

nt 
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<460  IP  (Z  <  2,5)  THEN  COTO  4500 

44?0  ERF:1! 

4480  GOTO  4660 

4490  'initalize  variables 

4500  TEST:1 

4510  BRF:0! 

4520  ANS:0! 

4530  NliO 

4540  ACT:1 

4550  IP  (NUO)  THEN  GOTO  4600 

4560  INC:N1 

4570  ACT:ACTtINC 

4580  INCzINC-l 

4590  IP  (INCOO)  THEN  GOTO  45Î0 

4600  EBS:(((-l!r(Nl))»((Zr(2!»MUl!)))/((2!«flltl!)»ACT) 

4610  ANSrRESUNS 

4620  NUNltl 

4630  IF  (8BS>:(TBST/1E+11))  THEN  GOTO  4540 

4640  IP  ((EBSt(-lM>:(TBST/lB+ll))  THBH  GOTO  4540 

4650  ERF:ANS»2/(3.1415927r(.5)) 

4660  IF  (NBGiO)  THEM  GOTO  4680 

4670  BEP:ERPt(-l) 

4680  'RINT  X.BBP 

4690  RETURN 

4700  'subroutine  for  screen  display 

4710  PBiNT°t>tn>>ttttttttntututttttstttnttnustnntnttttntttttttnttnt 

mtt" 

4720  PRINT  "tttmtttttttnttnnttttttntSPILL  PROFILE«»ttt»t»lUJ»ut»ntunt 

jntj" 

4730  PBINT'tttttnnttttnutttttntttttntnittttttttttutntttnntstttttttttt 

ttm" 

4740  PRINT 

4750  'insert  decay  rate  if  applicable 

4760  PRINT  "SUBSTANCE  SPILLED  WAS:  "  CHEH$  ;IP  DFLAdl  THEN  PRÎHT  "DECAY  RATE  I 

S"  DBCAÏ  "S'-l" 

4770  PRINT  "IN  THE  "  RIVEB$ 

4780  PRINT  "THE  SOURCE  OP  THE  SPILL  «AS:  "  SOURCE) 

4790  PRINT 

4800  PRINT  'TIHB  OP  SPILL:  FROM  "  START)"  TO  "SPEND$ 

4810  PRINT 

4820  IP  FLAGrl  THEN  GOTO  4850 

4830  PRINT  "THIS  PROFILE  IS  AT  A  DISTANCE  OF  "DIST(l)"  leters  DOHNSTRBAH  FROM  TH 

E  SPILL 

4840  'alert  if  portion  of  spill  truncated 

4850  PRINT 

4860  IF  TRUNC:!  THEN  PRINT  "ONLY  THE  PORTION  OP  THE  SPILL  WHICH  HAS  OCCURED  BBFO 

RE  THE  TIHB  IN  WHICH  BACH   CONCENTRATION  HAS  OCCURBD  IS  USED" 

4870  PRINT 

4880  PRINT  "NUMBER", ■  TIME  FROM  START", "CONCENTRATION"  ; 

4890  'insert  distance  to  concentration  if  discrete  concentrations  are  used 

4900  IF  FLAG^l  THEN  PRINT  "    DISTANCE" 

4910  IF  PLAG:0  THEN  PRINT 

4920  PRINT  ■  OP  SPILL  (hrs)","  ag/l"; 

4930  IF  PLAGil  THEN  PRINT  "  i   " 

4940  IF  PLAG^O  THEN  PRINT 

4950  PRINT 

4960  FOR  1:1  TO  NJ 

4970  'goto  decay  subroutine  if  applicable 

4980  IP  DFLAG:1  THEN  GOSUB  5610 

4990  PRINT  USING  'mi',   I  ; 
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5000  PRINT  ■      •; 

5010  PEINT  USING  'HU\Hi.Ht',    TT(I1; 

5020  PEINT  "  '; 

5030  PRINT  USING  'MM ;  TCON(I); 

5040  IF  FLAGiO  THEN  PEINT 

5050  IF  FLAG:1  THEN  PEINT  "      "DIST(I) 

5060  NEXT  I 

5070  INPUT  "DO  YOU  ÏISH  TO  PRINT  THIS  TABLE  (Y/N)  ";ZX1$ 

5080  IP  im-.'V   ÛE  ZXI$:"y"  THEN  PEFLAGxl  :GÛSUB  5140 

5090  IF  UH-'r   OR  ZUi-.'n'   THEN  GOTO  5130 

5100  'go  back  to  option  iienu  after  printing 

5110  IP  PRFLAG^l  THEN  GOTO  5130 

5120  BEEP:GOT0  5070 

5130  PRFLAGrO  :EETUEN 

5140  'subroutine  for  printer 

5150  LPEINT''ntntntntttttttttttttntsttttt>sttntnttutnttttnntnt{ttutt 

UtJtt" 

5160  LPRINT  'tuttttnttnintuunmnnsPiLL  PROFILBtntJunn«m»ntnm 
ttnn" 

5170  LPBINT''<nt:nntmttntt tut tttn tsnu ttt nut ttntint nut tnnntnss 
tittn' 

5180  LPEINT 

5190  LP8INT  "SUBSTANCE  SPILLED  WAS:  '  CHEH$  :IF  DPLAG:1  THEN  LPRINT  "DECAY  RATE 

IS"  DECAY  "3*-r 

5200  LPRINT  "IN  THE  "  RIVER} 

5210  LPEINT  "THE  SOURCE  OF  THE  SPILL  «AS:  "  SOURCE} 

5220  LPRINT 

5230  LPRINT  "TIME  OF  SPILL:  FROM  "  START}"  TO  "SPEND} 

5240  LPRINT 

5250  IP  FLAG:1  THEN  GOTO  5290 

5260  LPRINT  "THIS  PROFILE  IS  AT  A  DISTANCE  OF  "DIST(l)"  aeters  DOKNSTREAH  FROH  T 
HE  SPILL 

5270  LPRINT 

5280  IF  TRUNC:1  THEN  LPRINT  "ONLY  THE  PORTION  OF  THE  SPILL  WHICH  HAS  OCCURED  BEF 

ORE  THE  TIME  IN  «HIGH  BACH   CONCENTRATION  HAS  OCCURED  IS  USED" 

5290  LPRINT 

5300  LPRINT  "NUMBER","  TIME  FROM  START", "CONCENTRATION"  ; 

5310  IF  FLAGrl  THEN  LPRINT  "    DISTANCE" 

5320  IF  FLAG:0  THEN  LPRINT 

5330  LPRINT  "  OF  SPILL  (hrs)","   nig/1"; 

5340  IF  FL.AG:1  THEN  LPRINT  "  m   " 

5350  IF  FLAGrO  THEN  LPRINT 

5360  LPEINT 

5370  FOR  1:1  TO  NJ 

5330  LPRINT  USING  "»»»";  I  ; 

5390  LPRINT  "      '; 

5400  LPRINT  USING  "Himti.m",    TT(I); 

5410  LPEINT  "  "; 

5420  LPEINT  USING  "MM ;  TCON(I); 

5430  IF  FLAGrO  THEN  LPRINT 

5440  IF  FLAG:1  THEN  LPRINT  "      "DIST(I) 

5450  NBXT  I 

5460  RETURN 

5470  'subroutine  to  calculate  inconplete  spills 

5480  'set  flag 

5490  TRUNCM 
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5500  TF(0):0  :THASS(0):0 

5510  'keep  original  tine 

5520  TBHP5:TP(I) 

5530  N:r 

5540  'set  partial  tine 

5550  TF(N)iT(Il) 

5560  'set  tiae  lB-9  hours  larger  than  spill  time 

5570  T(Il):T(Il)+lB-09 

5580  'linear  truncation  of  last  partial  spill  fluï 

5590  mSS((():TNASS(M)t(TF(N)-TF{[-l))/(TEHP5-TF(r- 
5600  RETURN 

5610  'DBCAY  SATE  SUBEOUTINB 

5620  TCOM(I):TCON(I)lBXP(-ltDBCAYtTT(I)J3600) 
5630  BBTUBN 
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